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ABSTRACT 


Identifying the influence of eastern boundary current (EBC) dynamical processes 
on sound speed structure can provide better ocean models for acoustic predictions in 
littoral regions. This study will explore the effects of currents, wind and eddies on the 
sound speed structure of two different EBC models, the North Canary Current System 
(NCCS) and the Leeuwin Current System (LCS). These systems represent classical 
features of EBCs as well as regional anomalies. This study seeks to introduce sound 
speed analysis to the sigma coordinate primitive equation models and determine which 
regions of EBCs experience the largest changes in sound speed and most intense 
gradients. Results of model runs show that the dynamics of EBCs lead to large changes 
in sound speed and distort the vertical sound speed profde. The greatest change in sound 
speed in either region is caused by upwelling in the NCCS. Surface and associated 
subsurface eddies in the LCS are the largest scale feature in the study. The undercurrent 
of the NCCS and Meddies present the most intense (horizontal and vertical) gradients of 
sound speed change. 
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I. INTRODUCTION 


As Navy operations continue to move into coastal regions, the need for accurate 
predictions of shallow water acoustics is of particular importance. Sound speed structure 
in “blue water” has been studied for many years and reasonable models exist to predict 
large scale and small scale sound speed variability. However, the complex ocean 
dynamics of littoral regions complicate the study of acoustics near the coast. Identifying 
the influence of dynamical processes like upwelling, shelf eddies, filaments, meanders, 
etc., on sound speed structure can provide better ocean models for acoustic predictions in 
littoral regions. 

This study will explore the effects of currents, wind and eddies on the sound 
speed structure of two different eastern boundary current (EBC) regions. The EBCs of 
specific interest are the North Canary Current System (NCCS) off Morocco and the 
Iberian Peninsula, and the Leeuwin Current System (LCS) off Western Australia. These 
systems represent classical features of EBCs as well as regional anomalies. Major 
subtropical EBCs such as the California, Peru, and Benguela Currents are wind-driven 
and characterized by climatologically weak (<10 cm/s), broad (-1000 km wide), cold 
surface flow towards the equator. Surface currents flow in the direction of the prevailing 
winds with a poleward undercurrent below. The EBC thermocline is typically shallow 
(<30 m depth) and the region has high biological production due to vast regional 
upwelling (Parrish et ah, 1983). While the NCCS and LCS contain a number of these 
classical features, both regions include notable anomalies. 


A. NORTH CANARY CURRENT SYSTEM 

The NCCS stretches from ~10°N to ~45°N along the coast of northwestern Africa 
and the Iberian Peninsula. This system is considered a classical EBC and marks the 
closing eastern boundary of the North Atlantic Gyre. Typical of other EBCs, the mean 
equartorward Canary Current (CC) is a relatively broad (-1000 km), slow (-10-30 cm/s), 
yearlong surface flow extending from depths of -500 m (Wooster and Reid, 1976). A 
poleward undercurrent exists near the coast beneath the CC and is relatively narrow 


1 



(-10-40 km) and weak (-“2-10 cm/s). The undercurrent is strongest between depths of 
-100 m and -600 m, varying with season and latitude. In winter, the CC shoals to the 
north near Cabo da Roca and fonns a component commonly referred to as the Iberian 
Current (IC) (Haynes and Barton, 1990). The IC is a relatively narrow (-25-40 km), 
weak (-20-30 cm/s), seasonal current found trapped near the coast against the shelf break 
(Fiuza, 1980; Frouin et al., 1990; Haynes and Barton 1990). 

The NCCS is predominately influenced by equatorward, upwelling favorable 
winds produced by the Azores High. The Azores High is a semi-permanent subtropical 
high pressure system located in the northeastern Atlantic Ocean, similar in nature and 
behavior to the North Pacific Subtropical High (Nelson, 1977). The center of the high 
migrates meridionally with the seasons, causing wind stress values over the NCCS to 
vary temporally. 

Observations of the sea surface in the NCCS region have shown highly energetic 
mesoscale features, including jet-like surface currents, meanders, eddies and filaments. 
Satellite sea surface images have shown nearshore upwelling during periods of upwelling 
favorable winds with several narrow filaments of cooler water extending off the coast of 
the Iberian Peninsula (Fiuza and Sousa, 1989) and Cape Ghir in northwest Africa (Van 
Camp et al., 1991; Hagen et al, 1996). In these images, upwelling filaments often extend 
-80-150 km offshore with alongshore spacing of -80-100 km between filaments (Haynes 
et al., 1993), and tenninate with dipole eddy pairs (Fiuza et al., 1982; Barton, 1998). 
Observations have also shown anti-cyclonic and cyclonic pairs of mesoscale eddies on 
the order of 100 km off the coast of the Iberian Peninsula (Fiuza, 1984; Stammer et al., 
1991). These mesoscale features have been observed during periods of upwelling- 
favorable winds and are located near coastline irregularities such as capes. These 
observations provide evidence that wind forcing along with coastline irregularities appear 
to be important mechanisms in the fonnation and continuance of many of the mesoscale 
features found in the NCCS domain and other EBC regions (Batteen et al., 2000). 

A unique feature that distinguishes the NCCS from other classical EBCs is the 
influence of the Mediterranean Outflow (MO). Due to its high density, the MO flows 
down the continental slope in the Gulf of Cadiz as a narrow gravity current and entrains 
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the less dense North Atlantic Central Water. Coriolis acceleration causes the outflow to 
deflect to the right as it leaves the strait (Bower et al., 1997). The highly saline plume of 
the MO dilutes, thickens, and becomes vertically differentiated into two distinct cores as 
it follows the northern rim of the Gulf of Cadiz (Iorga and Lozier, 1999; Bower et al., 
1997). At ~7°W in the Gulf of Cadiz, the two cores exist at depths of -600-900 m and - 
1100-1200 m. Both cores continue to flow westward along the southern coast of Spain 
and turn poleward around Cabo de Sao Vincente (Amber and Howe, 1979; Iorga and 
Lozier, 1999). A third, shallower, poleward core of Mediterranean water has also been 
traced from the Strait of Gibraltar northward to ~38.5°N off western Portugal (Amber, 
1983). Additionally, climatological cyclonic circulation in the southwestern Gulf of 
Cadiz spreads the salty MO south ~34°N (Iorga and Lozier, 1999). 

In the past three decades, studies on how the Mediterranean salt tongue is 
maintained have prompted a search within the region for submesoscale coherent vortices 
(SVCs) or Meddies (e.g., Armi et al., 1989). Observations reveal that Meddies are 
lenslike features, containing a core of Mediterranean waters 200-1000 m thick centered at 
about 1000 m, which are less stratified than the surrounding fluid. They are typically 20- 
100 km in diameter and have temperature and salinity anomalies of approximately 2°C 
and 1 psu. This density structure supports an energetic anti-cyclonic rotation with peak 
azimuthal velocities of about 30 cm/s and rotation periods of 2-7 days (Armi et al., 1989; 
Hebert et al., 1990). The primary generation region of Meddies is widely accepted to be 
near Cabo de Sao Vincente off southwest Portugal. Several different trajectories of 
Meddies have been observed, including a southwestward movement into the Canary 
Basin, and westward translations south of the Azores (Richardson and Tychensky, 1998). 

B. LEEUWIN CURRENT SYSTEM 

The Leeuwin Current system is a sharp contrast to the typical EBC equatorward 
surface currents and poleward undercurrents of the NCCS. The Leeuwin Current is an 
anomalous, surface EBC off the coast of Western Australia that flows poleward along the 
continental shelf (~200 m depth) from Shark Bay to Cape Leeuwin. At the cape, it 
swings eastward as far as the Great Australian Bight (Cresswell and Golding, 1980). 
Like the NCCS region, prevailing winds along the coast of Western Australia are 
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predominantly equatorward. Unlike the NCCS, the Leeuwin Current is driven by an 
uncharacteristically large thennohaline gradient that overwhelms the winds stress. There 
is general agreement (Godfrey and Ridgway, 1985) that the LCS is generated by a 
meridional pressure gradient resulting from excessive heating in the equatorial region and 
large amounts of cooling in the poleward region. The source for the Leeuwin Current is 
predominantly an alongshore steric height gradient due to tropical Pacific water from the 
Indonesian throughflow (Godfrey and Ridgway, 1985; Hirst and Godfrey, 1993) and is 
augmented by a source from the North West Shelf (e.g., Gentilli, 1972), possibly having 
its origin in the Pacific Ocean (Hirst and Godfrey, 1993). As a result of this strong 
inflow of subtropical waters towards the coast, the LCS intensifies poleward (Batteen et. 
al„ 1992). 

Observations off the coast of Western Australia have shown the LCS is a strong 
(>150 cm/s at times), narrow (<100km wide), poleward surface current that opposes the 
prevailing wind direction (Godfrey and Ridgway, 1985; Cresswell and Golding, 1980; 
Godfrey et ah, 1986). Associated with the system are near-surface isopycnals that slope 
downward, surface dynamic height fields that increase toward the coast, and an 
equatorward undercurrent (Thompson, 1984; McCreary et al., 1986). Additionally, the 
LCS contains anomalous warm water at the surface and a deep (>50 m depth) 
thennocline (Thompson, 1984). The warm flow is credited with the appearance of 
tropical marine species off both the western and southwestern Australian coasts (Pearce 
and Griffiths, 1991) and lower biological production due to vast regions of downwelling 
(Batteen et. al., 1992). 

Below the Leeuwin Current is an anomalous equatorward undercurrent (e.g., 
Church et al., 1998). Current meter data from the Leeuwin Current Interdisciplinary 
Experiment (LUCIE) showed the undercurrent to be narrow and confined to the 
continental slope between ~250 m and 450 m depth (Smith et al., 1991). Between -250 
m and 350 m depth, the undercurrent can attain speeds comparable to the surface flow, 
exceeding on average 10 cm/s, up to ~30 cm/s (Smith et al. 1991). Although there is 
evidence of a 20 cm/s westward flow centered between ~400 m and 700 m depth off the 
southern coast (~l I6°E), no studies to date have clearly resolved whether there is a 

westward undercurrent off Southern Australia (Cresswell and Peterson, 1993). 

4 



From February to May each year, the Leeuwin Current strengthens considerably 
and strong eddies develop. Other dynamical features that have been observed in the LCS 
are offshoots of the current, meanders and both anticyclonic and cyclonic eddies (e.g., 
Legeckis and Cresswell, 1981). 


C. EFFECTS ON SOUND SPEED 

Sound speed profdes may be varied by a number of local conditions. Features 
such as internal waves, eddies or tides create variations in water density with depth 
(vertical displacement of the ocean layers). Near polar regions or river estuaries, profiles 
are affected by horizontal advection of water masses (Lurton, 2002). The classical and 
unique features of the NCCS and LCS change sound speed through both processes, 
vertical displacement or horizontal advection. Surface and subsurface eddies in both 
regions impact sound speed by their vertical displacements. Upwelling is another form of 
vertical displacement where deep water masses are lifted to the upper ocean. The 
anomalous surface current of the LCS advects wann water along the coast, causing 
intense horizontal gradients of sound speed. Mediterranean waters exiting the Strait of 
Gibraltar into the North-East Atlantic is a classic example or water mass exchange and 
horizontal advection (Lurton, 2002). Through these processes, the features of both 
regions will affect sound speed. 

Previous numerical studies by Batteen et al. (2006) investigated the effects of 
annual climatological wind forcing, initialized thennohaline gradients, and topography on 
the Leeuwin Current system. Similar studies have been conducted for the NCCS region, 
investigating the roles of wind forcing, bottom topography and thermohaline gradients on 
classical as well as unique features of the Northern Canary Current system (Martinho, 
2001). This study seeks to extend prior efforts in these areas by introducing sound speed 
analysis. The analysis will determine which regions of EBCs experience the largest 
changes in sound speed and most intense gradients. A comparison of the NCCS and LCS 
region will provide a diverse compilation of classical and unique features and their 
effects on sound speed. Knowledge gained from this study will allow Navy 
oceanographers to better understand the complex ocean dynamics observed in EBCs and 
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more accurately predict sound speed variability in these regions. In turn, the Undersea 
Warfare (USW) community will be better equipped to support submarine operations in 
these areas. 
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II. MODEL DESCRIPTION 


A. DATA SETS 

The model was initialized using annual climatological temperature and salinity 
maps obtained from Levitus and Boyer (1994) and Levitus et al. (1994). These data sets 
use a 1 by 1 degree horizontal resolution at thirty-three vertical levels shown in Table 1. 

Topographic data was obtained from the Institute of Geophysics and Planetary 
Physics, University of California San Diego (Sandwell and Smith, 1996). The bottom 
topography has a resolution of 2 minutes (1/30 of a degree) and was complied from 30 
years of bottom echo soundings by ships. Gravity data was used to estimate soundings in 
the vast regions in which there was sparse or no ship data. 

Climatological wind fields were obtained from the European Centre for Medium 
Range Weather Forecasts (ECMWF) near-surface wind analyses (Trenberth et al., 1990). 
The wind data uses a 2.5 by 2.5 degree grid. 


B. PRE-PROCESSING 


The original topography was interpolated to the resolution of the POM model for 
each current system. In the NCCS model, the resolution was 3 by 3.7 km in the shallow 
regions and 6 by 7.4 km in deeper region, resulting in a total 287 by 241 horizontal grid. 
The variable grid in this model allows the use of more (less) points in regions of large 
(small) gradients. The highest resolution was used where the values of the slope 


parameter (defined by Mellor et al., 1998, as 


1 SH\ 

2 *H 


where H is the average depth and 


SH is the difference in depth between two adjacent cells) were the largest in both the 
latitude and longitude directions. The maximum suggested slope parameter to be used in 
sigma coordinate models is 0.2 (Mellor, 1996). Since the slope parameter was greater 
than this maximum value over much of the area, the topography was smoothed with a 
linear, 2D low-pass filter in order to meet this criterion. The new depth of each point 
calculated with this filter was a non-weighted average of 15 by 15 points surrounding the 
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point. Depths greater than 2500 m were reassigned to the maximum depth of 2500 m. 
Land was assigned a depth of 10 m to avoid division by zero in the model and the Strait 
of Gibraltar was closed (Figure 1). 

Similarly, the LCS model horizontal grid resolution varies from 5.7 km by 6 km 
to 13.7 km by 10.2 km, resulting in 298 by 233 points. The finest resolution, 5.7 km by 6 
km, was located over the continental shelf and slope around the southwest corner of 
Australia; the bottom slope and temperature gradient were greatest in this area. Original 
topography was smoothed with a one dimensional (ID), iterative method (Martinho and 
Batteen, 2006). This filtering process has been shown to satisfy the slope criterion while 
maintaining coastline irregularities, continental shelves, and relative maxima such as 
seamounts and islands. To avoid altering slope parameters after they have been adjusted, 
a two dimensional (2D) averaging filter was first applied to the coastline. For the 
topography, the signed slope parameter was then calculated along each grid line in one, 
initial direction over the domain. Any two adjacent cells with a slope parameter greater 
than the limit were adjusted to the limit value, and the line was recalculated. When all 
grid lines in the initial direction were adjusted, the topography matrix was rotated 90 
degrees. Each process was repeated until the topography was adjusted for all directions 
(rotated 360°). To avoid division by zero in the model, land was assigned a depth of 20 
m and depths greater than 2500m were reassigned to 2500 m (Figure 2). 

In both models, the annual climatological temperature and salinity values were 
interpolated for the horizontal spatial resolution of the model and for the 21 vertical 
sigma levels with a three dimensional (3D) linear interpolation scheme. Daily seasonal 
winds were averaged over time in order to obtain the annual, non-weighted, average wind 
vector field. The wind vectors were also interpolated for the horizontal spatial resolution 
of the model with a two dimensional linear interpolation scheme. Since wind stress 
rather than the wind forces the model, the annual wind vectors were then used to compute 
the wind stress. 
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C. BRIEF MODEL DESCRIPTION 

The Princeton Ocean Model (POM) was used in the model studies of the NCOS 
and LCS. The POM, a well-documented model, is a primitive equation, free surface 
model with a second-moment turbulence closure scheme (Mellor and Yamada, 1982). 
Using bottom-following sigma levels, the POM can realistically simulate processes 
associated with continental shelf flows and bottom boundary layer dynamics in local 
domains such as bays, estuaries and coastal regions. This model utilizes 21 vertical 
sigma levels. The sigma values ranged from zero at the surface to minus one at the 
bottom with the vertical grid spacing proportional to the ocean depth. The vertical 
resolution was higher near the surface and the bottom in order to resolve the boundary 
layers in these regions. An implicit vertical time differencing scheme was used to 
eliminate the time constraint of the vertical grid near the surface, bottom and in shallow 
water. The sigma values for both models are found in Table 2. 

The prognostic variables of the model are potential temperature, salinity, density, 
three velocity components, surface elevation, turbulent kinetic energy and length scale. 
The model has a split time step interval for external and internal modes. The external 
mode solves the equations for the vertically integrated momentum equations. It also 
provides the sea surface and barotropic velocity components and has a time step of six 
seconds. The internal mode solves the complete 3D equations and has a time step of 300 
seconds in the LCS model and 180 seconds in the NCCS model. 

A Smagorinsky formulation (Smagorinsky et ah, 1965) was used for the 
horizontal diffusion in which the horizontal viscosity coefficients depend on the grid size, 
the velocity gradients and a coefficient. Consistent with other POM studies, a value of 
0.2 was assigned to this coefficient (Ezer and Mellor, 1997). 

D. INITIALIZATION, FORCING AND BOUNDARY CONDITIONS 

The model was initialized with annual temperature and salinity values obtained 
from Levitus and Boyer (1994) and Levitus et al. (1994). Since the model runs reached a 
quasi-equilibrium state in a relatively short time (<60 days), zero salinity and temperature 
fluxes have been prescribed at the ocean surface. 
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The model was forced from rest with the annual ECMWF wind fields which were 
interpolated for the model grid. In the NCCS region, the wind stress is stronger in the 
southern region of the domain and weaker off Iberia and in the Gulf of Cadiz (Figure 3). 
In the FCS region, westerly winds dominate the southern portion of the domain (Figure 
4). Along the west coast of Australia the winds are southerly to south-southeasterly and 
generally increase in strength away from the coast. 

Correct specification of the open boundary conditions (BCs) is very important to 
achieving realistic results without reflections, clamping, spurious currents or numerical 
alteration of the total volume of water in the model. Generally, boundary conditions are 
specific to each model and study. For models with a free surface, such as the two used 
here, BCs should be transparent to the waves. A gradient boundary condition (Chapman, 
1985), which allows geostrophic flow normal to the boundary, works best for the 
elevation. For baroclinic velocity components nonnal to the boundary, an explicit wave 
radiation scheme based on the Sommerfeld radiation condition was used. For inflow 
situations, the model was forced with annual temperature and salinity values obtained 
from Fevitus and Boyer (1994) and Fevitus et al. (1994), while an advection scheme was 
used in outflow situations. An improvement on previous models is the introduction of an 
advanced volume constraint subroutine based on research perfonned by Marchesiello et 
al. (2001), which drastically reduced the loss of volume at the boundaries to insignificant 
levels. 

For the barotropic velocity components, a Flather radiation plus Roed local 
solution (FRO) was used. Palma and Matano (2000) showed good results with the FRO 
during BC tests to determine the BCs response to an alongshelf wind stress. Palma and 
Matano (1998) also showed that the FRO BC had good reflection properties and results in 
a test that determined the BC response to the combined action of wind forcing and wave 
radiation. These tests were executed with the barotropic version of POM and compared 
with benchmark results (no boundary conditions). 
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E. SOUND SPEED CALCULATION 

Sound speed in this study was calculated using the Chen and Millero equation in 
the UNESCO algorithm (Fofonoff and Millard, 1983), and is a function of temperature, 
salinity, and pressure (depth and latitude). Equations for the speed of sound in seawater 
based on laboratory measurements have been published by Wilson (1960), Del Grosso 
(1974) and Chen and Millero (1977). The equation of Chen and Millero is generally 
accepted as the sound speed standard, although there is much debate in scientific 
literature about the accuracy and range of applicability of these equations. Studies have 
shown the difference between Del Grosso and Chen and Millero sound speed calculations 
increases with depth to almost 0.6 m/s at 5000 dbar for realistic ocean temperatures and 
salinities (Dushaw et ah, 1993). Since the models in this study extend to 2500 m depth, 
the Chen and Millero equation was selected. Depth was converted to pressure by the 
Saunders (1981) depth-to-pressure conversion technique. 
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III. RESULTS 


A. NORTH CANARY CURRENT SYSTEM MODEL 

Initial climatology of the NCCS region can be seen in Figure 5 and Figure 6. 
Temperature and salinity decrease roughly linearly from north to south by 3°C and 0.8 
psu, respectively. Surface sound speed (Figure 7) increases by -5 m/s from north to 
south and is relatively homogeneous, with the exception of waters north of 37°N which 
resemble the temperature gradient. A cross-sectional profile of salinity at day 1 through 
the gulf at 36°N (Figure 8) shows the Mediterranean Outflow (MO) signature at -1200 m 
depth with values of -36.2 psu, and characteristically high salinity values of the North 
Atlantic waters in the upper 300 m. Evidence of the outflow waters can be seen across 
the entire region (not shown), although the highest values of salinity are found in the gulf. 
Cross-sections of salinity at 32.5°N (Figure 9) and 39.3°N (Figure 10) show that values at 
depths of the MO decrease by -0.4 psu to the south and -0.1 psu to the north of the gulf 
(compare Figures 8, 9, and 10). The MO signature appears in sound speed profdes for 
the same day and latitudes (Figure 11), creating two distinct sound channels at depths of 
approximately 500 m and 2000 m. This particular profile is unique to this region where 
the warm, salty waters of the Mediterranean mix with the North Atlantic. Important 
acoustically, the double minimum features are strongest in the gulf and decrease towards 
the north and south, exemplifying the same spreading of the Mediterranean waters seen in 
salinity cross sections (Figures 11, 12 and 13). 

1. Upwelling 

Using these initial conditions and the winds, the model was run for 90 days. By 
day 20, upwelling is evident along the coast (Figure 14). Cross sections profding the 
difference in sound speed from the initial conditions indicate a decrease over time in 
sound speed caused by the cooler waters. For example, the effect on sound structure at 
latitude 39.3°N (Figure 15) has decreased -15 m/s. The upwelling extent continues to 
broaden along the entire coastline through day 45 and reaches its maximum extension by 
day 63 (Figure 16). This wide and shallow feature subsequently decreases sound speed 
along the shelf by a maximum of -19 m/s at all latitudes from the surface to -200 m 
depth in the north (Figure 17), -400 m depth in the gulf (Figure 18), and -500 m depth in 
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the south (Figure 19). Although upwelling primarily impacts the mixed layer, the 
horizontal gradient of sound speed in the upper layer is not large and is changing -19 m/s 
over -280 km. 

2. Undercurrent 

A warm, poleward undercurrent appears by day 20 north of the Gulf of Cadiz. 
The undercurrent reaches its maximum intensity by day 36, after which eddies and 
Meddies begin to fonn and spin off. South of the gulf, the sound speed increases by 5 
m/s due to the warm waters of the undercurrent flowing northward (Figure 20). 
Similarly, cross sections within the gulf (36°N) also show a 5 m/s increase (Figure 21). 
North of this latitude, the undercurrent advects MO waters and has a greater effect on 
sound speed structure. A cross section of sound speed difference at 39.3°N shows a 9 
m/s increase due to the added salinity of the MO to the wann undercurrent (Figure 22). 
These changes in sound speed occur over a shorter longitudinal distance than upwelling, 
extending -15 km from the point of maximum increase to the point of no change in 
sound speed. The resulting horizontal gradient of sound at depth along the slope is 
significantly larger than in the mixed layer. 

3. Meddies 

In the undercurrent around Cabo Roca (39.3°N), a Meddy forms by day 45, 
generating a dipole pair (Figure 23). Sound speed across the Meddy increases by -6 m/s 
within 50 km horizontally and 600 m vertically (Figure 24). As the Meddy mixes MO 
waters with the North Atlantic, its effect on sound speed weakens and the gradients 
become less condensed. By day 63, the Meddy and dipole pair are still visible in the 
39.3°N cross section, but show only a 3.5 m/s increase over 50 km (Figure 17). To the 
south, an eddy appears in the undercurrent at 32°N by day 63; however, this eddy forms 
after the undercurrent has been established and the MO waters have already been 
advected to the north. Consequently, this eddy does not contain the MO signature and so 
is not considered a Meddy. Without the added salinity of the MO, the sound speed 
increases only -2 m/s across the -100 km feature (Figure 25). 

4. Filaments 

On the surface, the NCCS generates several filaments of cooler water from 
upwelling. By day 63, a topographic feature at 33.9°N causes a surface eddy to spin off 
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from these filaments. The well-defined, cold-core eddy causes a decrease in sound speed 
of -12.5 m/s across -40 km and reaches below the mixed layer to a depth of 1500 m 
(Figure 26). Similar to the undercurrent eddy at 32°N, this eddy is not influenced by the 
MO, despite affecting depths deeper than the outflow. 

5. Classical and Unique Features 

The NCCS was chosen for this study for its classical upwelling and equatorward 
surface current as well as the unique impact of the Mediterranean Outflow. Table 3 
summarizes the effects on sound speed by each feature in this region. Results are 
interpreted in terms of effective change in sound speed from initial conditions as well as 
spatial extent. Of critical acoustical importance are the spatial gradients of sound speed 
caused by these features. 

In the NCCS, upwelling clearly incurs the greatest change on sound speed, 
affecting a large surface area through the mixed layer. Sound speed fluctuation in this 
feature is driven by the change in temperature and reaches a maximum at the coast soon 
after the appearance of upwelling, remaining constant while upwelling extent continues 
to grow steadily. As expected, upwelled waters are significantly cooler than initial 
surface conditions and cause a decrease in sound speed. Consequently, eddies which spin 
off upwelling filaments also cause a large decrease in sound speed. These mesoscale 
surface eddies are significant acoustical features due to their large gradients of sound 
speed. Although the surface eddy affects sound speed -6.5 m/s less than upwelling, 
sound speed structure in the eddy changes over a much shorter longitudinal distance, 
leading to larger gradients. 

Effects of Mediterranean waters in the NCCS model can be seen in results of the 

northward undercurrent. Within the Gulf of Cadiz and south of the gulf, the undercurrent 

incurs a consistent change, primarily due to the increase in temperature from southern 

waters. Although the undercurrent narrows to the north, its influence on sound speed 

increases due to the added salinity of the MO advected northward in the current. Similar 

results are seen in a comparison of the undercurrent eddy and Meddy. The eddy found 

near the south boundary of the domain causes very little change to the sound speed 

structure. The size of the Meddy off Cabo da Roca is on the same order as the 

undercurrent eddy; however, the Meddy contains the MO signature. Consequently, the 
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Meddy incurs a greater sound speed change. The NCCS model shows the Meddy 
dissipating over time; however, Meddies in the North Atlantic play a fundamental role in 
maintaining the large-scale salt tongue and gradual decrease in salinity is expected 
(Bower et ah, 1997). Whereas temperature is a foremost factor affecting sound speed, the 
NCCS represents a region uniquely affected by salinity. 


B. LEEUWIN CURRENT SYSTEM MODEL 

Figures 27, 28, and 29 show the initial climatology for the LCS model. The 
meridional temperature gradient of this region (Figure 28) is notably larger than the 
NCCS region with a maximum temperature near Shark Bay of —24°C and a minimum of 
14°C at the southern edge of the model domain. The climatological surface salinity 
(Figure 27) shows a maximum of -38.9 psu at the center of the western edge of the 
model domain and a minimum of —35.3 psu along the northern Australian coast. Despite 
a slight increase in salinity between 34°S and 30°S (Figure 27), the surface sound speed 
decreases roughly uniformly from north to south. The resulting difference between 
maximum and minimum surface sound speed is -8 m/s (Figure 29). Cross-sectional 
profiles of sound at various latitudes (e.g., 28°S shown in Figure 30) indicate a minimum 
of-1490 to 1494 m/s over the entire region at -1200 m depth. 

1. Leeuwin (Surface) Current 

By day 20, the narrow (-90 km) Leeuwin Current has fonned primarily as a result 
of the thermal gradient over the region (Figure 31). Warm water is advected poleward 
within the narrow flow, resulting in large temperature gradients around Cape Leeuwin. 
Figure 32 shows large horizontal and vertical gradients across this feature with an 
increase of -9 m/s. The current quickly strengthens over time and, by day 36, has 
reached a width of-160 km notably without the formation of eddies (Figure 33). Sound 
speed increases by -12 m/s (Figure 34). After day 36, the current reduces in intensity 
and eddy generation begins. 

2. Undercurrent 

Similar to the surface current, an equatorward undercurrent appears around Cape 
Leeuwin by day 20 (Figure 35). By day 36, the undercurrent appears at all latitudes on 
the west coast (not shown) and has reached its maximum intensity without the formation 
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of eddies. Unlike the wide, offshore surface current, the undercurrent narrows 
considerably from -90 km horizontally and 1000 m vertically to -25 km and 500 m as it 
flows equatorward. Topography traps the current to the shelf at -1200 m depth. 
Consequently, this feature primarily affects the shelf break region, resulting in a -3 m/s 
to 7 m/s decrease in sound speed. 

3. Upwelling 

Although the surface current flows opposite to the wind direction, southerly winds 
generate minimal upwelling predominantly in the bight around 33.5°S and between 28°S 
and 30°S (Figure 36). Note that upwelling does not appear where the LCS flows close to 
the coast. In the regions that are affected by upwelling, water relatively cooler than the 
surface current appears by day 20; however, a cross section of the bight shows almost no 
change from the initial sound structure on the shelf (Figure 32). Due to the shallow depth 
of the shelf and the wann water advected by the LCS from the north, upwelled waters 
cause either no change or a slight increase from sound speed initial conditions. 
Upwelling reaches a maximum extent at day 63 and increases sound speed by 1 m/s in 
the north (28°) and 4 m/s in the bight. The waters are confined to the shelf, extending 40 
km to 70 km from the coast and from the surface to 100 m depth. 

4. Eddies and Filaments 

Strong eddies are among the most notable features of this model run. By day 36, 
baroclinic instability generated by increasing vertical shear causes the current to meander 
(Figure 33). By day 45, meanders at 28°S and 35.5°S become wann-core eddies 
propagating westward (Figure 37). Cross sections at 28°S (Figures 38 and 39) show that 
the eddies grow to -250 km in width from day 45 to day 63 and affects the upper 1500 m 
of the water column by -7 m/s. Note that no mesoscale features in the undercurrent form 
at this latitude. At 35.5°S, smaller (-130 km) eddies spin off of Cape Leeuwin with 
warm-core eddies at the surface and cold-core eddies in the undercurrent (Figure 40). 
Unlike the surface eddy at 28°S, this warm-over-cold core eddy affects the entire water 
column, increasing the upper 1200 m by a maximum -7 m/s and decreasing the lower 
1300 m by 5 m/s. Equatorward winds around the cape help to develop eddy dipole pairs 
which propagate westward. The third and strongest eddy fonns by day 54 at 33.5°S, also 
as a result of opposing winds and topography (Figure 41). Unlike eddies off the southern 
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coast of Australia, this -180 km feature does not propagate westward and only affects the 
water column to -400 m depth. The warm core is more defined than previously formed 
eddies and produces a strong sound speed gradient across -90 km (Figure 42). This 
surface and undercurrent eddy causes a -10 m/s increase and -2 m/s decrease 
respectively in sound speed. 

5. Classical and Unique Features 

Table 4 summaries the results of the LCS model run and details the impact of 
each feature on sound speed. The strong meridional pressure gradient off the coast of 
Australia produces features unique to the LCS region, particularly the surface current. 
This feature incurs the greatest change in sound speed due to poleward advection of 
warm waters. Consequently, sound speed increases within the current and the greatest 
change from initial conditions occurs in the coolest waters around Cape Leeuwin. The 
surface current tapers towards the cape, creating a narrow, shallow feature with a large 
effect on sound speed. Similarly, the undercurrent advects cooler water and narrows as it 
flows northward, causing a region of decreasing sound speed between 1000 m and 1400 
m depth. Both the surface and undercurrent remain close to the topography. Eddies in 
the surface current and in the undercurrent tend to form off the same topographic 
features, causing changes in the sound speed of the entire water column. These eddies 
are significant influences on sound speed structure as they affect large, isolated volumes 
of the domain by as much as 10 m/s. 

In the LCS, forcing from the thermal gradient overwhelms wind stress. 
Upwelling is confined to the latitudes where the current has drifted westward away from 
the coast, allowing minimal upwelling on the shelf. Due to the strong advection of wann 
waters in the surface current, upwelling on the shelf is wanner than initial conditions 
though relatively cooler than the surface current. As a result, upwelling incurs a weak 
increase in sound speed and affects a relatively small area. 

C. NCCS AND LCS COMPARISON 

The objective of this study was to assess which dynamics of EBCs have the 
greatest impact on the sound speed structure. The NCCS and LCS were chosen for their 
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unique characteristics, i.e., the Mediterranean Outflow in the Gulf of Cadiz and the 
anomalous, poleward surface flow off the coast of Australia. As shown in Tables 3 and 
4, these characteristics produce a variety of features and results. 

The greatest change in sound speed in either region was caused by upwelling in 
the NCCS. Although the NCCS and the LCS are both influenced by equatorward, 
upwelling-favorable winds typical of EBC regions, the resultant upwelling differs 
dramatically between the two systems. In the NCCS, upwelling is present at all latitudes, 
has an extension on the order of 100 km, and causes a predictable decrease in sound 
speed due to cold waters. Contrary to this intuitive sound speed decrease in NCCS 
upwelling, the confined regions of upwelling in the LCS show an increase in sound speed 
due to the proximity of the surface current. The Leeuwin Current flows along the shelf 
break, restricting upwelling to the shelf. Consequently, the extent of NCCS upwelling is 
also significantly larger than the LCS and affects a larger volume of the domain. 

Despite differences in upwelling, both the NCCS and LCS surface currents 
generated eddies of great acoustical significance. Both currents experienced large 
changes in temperature; the NCCS flow contained cool waters from upwelling and the 
LCS advected warm waters from the thermohaline gradient. Eddies from either system 
contained cores of the surface current temperature. Although greater changes in sound 
were found in the NCCS due to the intensity and extent of upwelling, the LCS model 
simulated more surface eddies than the NCCS. The LCS surface eddies affected a greater 
surface area and volume, evidence of a stronger system. Eddies in both areas affected 
more than half the water column, extending to -1500 m depth. 

Although the LCS is a stronger current system, greater changes in sound speed by 
the undercurrent were found in the NCCS, particularly to the north of the gulf, due to the 
presence of the Mediterranean Outflow. The NCCS undercurrent is also narrower and 
deeper than the LCS, affecting a smaller volume but resulting in very large horizontal 
gradients of sound speed along the shelf slope. Eddies generated by the undercurrents in 
both regions are small in diameter and effect less than a 5 m/s change on sound speed. 
Meddies of the NCCS incur more change than undercurrent eddies of the LCS due to the 
added salinity of the MO advected by the NCCS undercurrent. 
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IV. SUMMARY 


The objective of this process-oriented study was to build on previous studies and 
determine the effects of eastern boundary current (EBC) dynamics on sound speed 
structure. The Princeton Ocean Model, a bottom-following sigma coordinate model, was 
used to explore the classical and unique features of two EBCs, the North Canary Current 
system (NCCS) and the Leeuwin Current system (LCS). The NCCS model simulated the 
upwelling, equatorward surface current and poleward undercurrent found in classic 
EBCs, as well as surface and undercurrent eddies and Mediterranean eddies (Meddies). 
The LCS model simulated the anomalous poleward surface current and equatorward 
undercurrent that is unique to the region, and the generation of energetic eddies. Using 
the equation by Chen and Millero (1977), sound speed was calculated for each model run. 

In the NCCS, upwelling appeared by day 20 and grew to a maximum extent by 
day 63. The wide and shallow feature decreased sound speed by ~19 m/s on the shelf at 
all latitudes from the surface to a maximum depth of 500 m. A single surface eddy 
formed in the south from an upwelling filament, creating a well-defined, cold-core 
feature of decreasing sound speed. A warm, poleward undercurrent formed and reached 
maximum intensity without eddies by day 36. The undercurrent increased sound by 5 
m/s in the gulf and south, and 9 m/s north of Cabo Roca. After day 36, the undercurrent 
generated two eddies. The eddy south of the gulf had a weak effect on sound, while the 
eddy in the north affected south three times more than in the south. These differences in 
sound speed in the north and south may be attributed to the presence of the 
Mediterranean Outflow. Water, salty waters of the Mediterranean Outflow are advected 
northwards in the undercurrent, causing stronger fluctuations in sound speed. 

In the LCS model, forcing from the thennal gradient overwhelmed wind stress, 
causing the strong, narrow Leeuwin Current to form by day 20. The current reached 
maximum intensity without eddies by day 36. Unlike classical EBCs, wann water was 
advected poleward within the flow, incurring a ~12 m/s increase in sound speed. An 
equatorward undercurrent formed simultaneously, narrowing considerably from south to 
north and advecting cold water from the polar region. The undercurrent had a smaller 
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effect on sound speed than the surface current and caused a decrease in sound. Southerly 
winds generated upwelling on the shelf in regions where the surface current had deviated 
from the coast. Although upwelled waters were relatively cooler than the surface current, 
upwelling caused a weak increase in sound from initial conditions. Energetic eddies 
formed by day 45 and continued to appear and grow through day 63. Eddies were warm- 
core on the surface and cold-core in the undercurrent, typically affecting a large volume 
of the domain and the entire water column. These features increased sound speed from 
the surface to -1500 m depth, and then decreased sound speed to -2500 m depth. 

The dynamics of EBCs lead to large changes in sound speed (5-15 m/s) and 
distort the vertical sound speed profile. The greatest change in sound speed from initial 
condition in either region was caused by upwelling in the NCCS. Surface and associated 
subsurface eddies in the LCS were the largest scale feature in the study. The 
undercurrent of the NCCS and Meddies presented the most intense (horizontal and 
vertical) gradients of sound speed change, incurring a significant change in sound across 
a relatively small feature. Although the LCS is a stronger current system, greater changes 
in sound speed by the undercurrent were found in the NCCS, particularly to the north of 
the gulf, due to the presence of the Mediterranean Outflow. Recommendations for 
further research include systematically investigating of the effects of each forcing 
mechanism on sound speed. Longer model runs would allow statistical analysis of 
individual features, determining the space and temporal correlation of upwelling, eddies, 
filaments, etc. Acoustic fields and ray traces may be computed with additional vertical 
layers in the model. 
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Figure 1. NCCS smoothed topography with depth in meters after applying a linear 
two-dimensional low-pass filter. (From: Sandwell and Smith (1996)) 
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Figure 2. LCS smoothed topography with depth in meters after applying a one¬ 
dimensional, iterative filter. Depth is in meters. 
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Figure 3. NCCS wind stress vector and magnitude (in color) in Pascals calculated 
from annual climatological ECMWF winds (From: Trenberth et al., 1990) 
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NCCS, Day 1, Surface. 
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Figure 5. NCCS, Surface temperature (°C) and velocity vectors (cm/s) on day 1. 
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NCCS, Day 1, Surface. 
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Figure 6. NCCS, Surface salinity (psu) and velocity vectors (cm/s) on day 1. 
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NCCS, Day 1, Surface. 
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Figure 7. NCCS, Surface sound speed (m/s) and velocity vectors (cm/s) on day 1. 
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NCCS, Salinity Cross-section at 36N, Dayl. 
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Figure 8. NCCS, Salinity (psu) cross-section at 36°N on day 1. 
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NCCS, Salinity Cross-section at 32.5N, Dayl. 

36.8 

36.6 

36.4 

36.2 
36 

35.8 

35.6 

35.4 

35.2 

-13 -12 -11 -10 -9 -8 

Longitude 

Figure 9. NCCS, Salinity (psu) cross-section at 32.5°N on day 1. 
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NCCS, Salinity Cross-section at 39.3N, Dayl. 
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Figure 10. NCCS, Salinity (psu) cross-section at 39.3°N on day 1. 
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NCCS, C Cross-section at 36N, Day 1. 
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Figure 11. NCCS, Sound speed (m/s) cross-section at 36°N on day 1. 
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NCCS, C Cross-section at 32.5N, Day 1. 
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Figure 12. NCCS, Sound speed (m/s) cross-section 32.5°N on day 1. 
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NCCS, C Cross-section at 39.3N, Day 1. 
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Figure 13. NCCS, Sound speed (m/s) cross-section 39.3°N on day 1. 
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NCCS, Day 20, Surface. 



-12 -10 -8 -6 
Longitude 

Figure 14. NCCS, Surface temperature (°C) and velocity vectors (cm/s) on day 20. 
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NCCS, Delta C Cross-section at 39.3N, Day20. 



Figure 15. 


NCCS, Difference in sound speed (m/s) from initial conditions, cross- 
section at 39.3°N on day 20. 
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NCCS, Day 63, Surface 
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Figure 16. NCCS, Surface temperature (°C) and velocity vectors (cm/s) on day 63. 
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NCCS, Delta C Cross-section at 39.3N, Day63. 
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Figure 17. 


NCCS, Difference in sound speed (m/s) from initial conditions, 
section at 39.3°N on day 63. 
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NCCS, Delta C Cross-section at 36N, Day63. 
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Figure 18. 


NCCS, Difference in sound speed (m/s) from initial conditions, cross- 
section at 36°N on day 63. 
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NCCS, Delta C Cross-section at 32.5N, Day63. 


Figure 19. 


NCCS, Difference in sound speed (m/s) from initial conditions, 
section at 32.5°N on day 63. 
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NCCS, Delta C Cross-section at 32.5N, Day36. 
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Figure 20. NCCS, Difference in sound speed (m/s) from initial conditions, cross- 
section at 32.5°N on day 36. 
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NCCS, Delta C Cross-section at 36N, Day36. 



Figure 21. NCCS, Difference in sound speed (m/s) from initial conditions, cross- 
section at 36°N on day 36. 
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NCCS, Delta C Cross-section at 39.3N, Day36. 
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Figure 22. NCCS, Difference in sound speed (m/s) from initial conditions, 
section at 39.3°N on day 36. 
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Figure 23. NCCS, Temperature (°C) and velocity (cm/s) at 1200 m depth on day 45. 
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NCCS, Delta C Cross-section at 39.3N, Day45. 
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Figure 24. NCCS, Difference in sound speed (m/s) from initial conditions, 
section at 39.3°N on day 45. 
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NCCS, Delta C Cross-section at 32N, Day63. 
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Figure 25. 


NCCS, Difference in sound speed (m/s) from initial conditions, cross- 
section at 32°N on day 63. 
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NCCS, Delta C Cross-section at 33.9N, Day63. 



Figure 26. 


NCCS, Difference in sound speed (m/s) from initial conditions, cross- 
section at 33.9°N on day 63. 
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Figure 27. LCS, Initial conditions of surface salinity (psu). 
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LCS, Day 0, Surface. 
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Figure 28. LCS, Initial conditions of surface temperature (°C). 
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LCS, Day 0, Surface. 



110 112 114 116 118 120 122 

Longitude 


Figure 29. LCS, Initial conditions of surface sound speed (m/s). 
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Figure 30. LCS, Sound speed (m/s) cross-section at 28°S for initial conditions. 
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LCS, Day 20, Surface. 
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Figure 31. LCS, Surface sound speed (m/s) and velocity vectors (cm/s) on day 20. 
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LCS, Delta C Cross-section at 33.5S, Day 20. 
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Figure 32. LCS, Difference in sound speed (m/s) from initial conditions, cross-section 
at 33.5°S on day 20. 
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LCS, Day 36, Surface. 
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Figure 33. LCS, Surface sound speed (m/s) and velocity vectors (cm/s) on day 36. 
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LCS, Delta C Cross-section at 33.5S, Day 36. 
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Figure 34. LCS, Difference in sound speed (m/s) from initial conditions, cross-section 
at 33.5°S on day 36. 
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LCS, Delta C Cross-section at 33.5S, Day 20. 
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Figure 35. LCS, Difference in sound speed (m/s) from initial conditions, cross-section 
at 33.5°S on day 20. 
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LCS, Day 20, Surface. 
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Figure 36. LCS, Surface temperature (°C) and velocity vectors (cm/s) on day 20. 
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LCS, Day 45, Surface. 
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Figure 37. LCS, Surface sound speed (m/s) and velocity vectors on day 45. 
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Figure 38. LCS, Difference in sound speed (m/s) from initial conditions, cross-section 
at 28°S on day 45. 
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LCS, Delta C Cross-section at 28S, Day 63. 
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Figure 39. LCS, Difference in sound speed (m/s) from initial conditions, cross-section 
at 28°S on day 63. 
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LCS, Delta C Cross-section at 35.5S, Day 45. 
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Figure 40. LCS, Difference in sound speed (m/s) from initial conditions, cross-section 
at 35.5°S on day 45. 
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LCS, Day 54, Surface. 


Figure 41. 
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LCS, Surface sound speed (m/s) on day 54. 
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LCS, Delta C Cross-section at 33.5S, Day 54. 
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Figure 42. LCS, Difference in sound speed (m/s) from initial conditions, cross-section 
at 33.5°S on day 54. 
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Table 1. Vertical levels and depths used by Levitus and Boyer (1994) and Levitus et 
al. (1994). 


Level 

Depth (m) 

Level 

Depth (m) 

Level 

Depth (m) 

1 

0 

12 

300 

23 

1400 

2 

10 

13 

400 

24 

1500 

3 

20 

14 

500 

25 

1750 

4 

20 

15 

600 

26 

2000 

5 

50 

16 

700 

27 

2500 

6 

75 

17 

800 

28 

3000 

7 

100 

18 

900 

29 

3500 

8 

125 

19 

1000 

30 

4000 

9 

150 

20 

1100 

31 

4500 

10 

200 

21 

1200 

32 

5000 

11 

250 

22 

1300 

33 

5500 


Table 2. Values of sigma levels for NCCS and LCS models. 


Level 

Sigma Value 

Level 

Sigma Value 

1 

0 

12 

-0.61538 

2 

-0.00961 

13 

-0.69231 

3 

-0.01923 

14 

-0.76923 

4 

-0.03846 

15 

-0.84615 

5 

-0.07692 

16 

-0.92308 

6 

-0.15385 

17 

-0.96154 

7 

-0.23077 

18 

-0.98077 

8 

-0.30769 

19 

-0.99038 

9 

-0.38462 

20 

-0.99519 

10 

-0.46154 

21 

-1.0 

11 

-0.53846 
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Table 3. Impact of features on sound in the NCCS. 



Latitude 

Ac 

Lx 

Lz 

A. Maximum upwelling extension 

39.3°N 

-19 m/s 

250 km 

200 m 

at day 63 

36°N 

-19 m/s 

270 km 

400 m 


32.5°N 

-19 m/s 

300 km 

500 m 



Ac 

Lx 

Core 



depth 

B. Undercurrent at day 36 

39.3°N 

+9 m/s 

15 km 

1600 m 


36°N 

+5 m/s 

20 km 

1800 m 


32.5°N 

+5 m/s 

25 km 

1500 m 


Day 

Ac 

Lx 

Lz 

C. Meddy at 39.3°N 

45 

+6 m/s 

50 km 

600 m 


63 

+3.5 m/s 

50 km 

600 m 



Ac 

Lx 


D. Eddy (undercurrent) at day 63 


+ 2 m/s 

100 km 


E. Eddy (surface) at day 63 


-12.5 m/s 

40 km 
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Table 4. Impact of features on sound in the LCS. 


A. Maximum upwelling extension 
at day 63 


B. Surface current at day 36 


C. Undercurrent at day 36 


D. Eddy (surface) 


E. Eddy (undercurrent) 


Latitude 

Ac 

Lx 

Lz 

28°S 

+ 1 m/s 

70 km 

50 m 

33.5°S 

+4 m/s 

40 km 

50 m 

Latitude 

Ac 

Lx 

Lz 

28°S 

+5 m/s 

200 km 

1000 m 

31.5°S 

+9 m/s 

120 km 

700 m 

33.5°S 

+ 12 m/s 

160 km 

500 m 

Latitude 

Ac 

Lx 

Lz 

28°S 

-3 m/s 

25 km 

500 m 

31.5°S 

-7 m/s 

25 km 

700 m 

33.5°S 

-4.5 m/s 

30 km 

800 m 

35.5°S 

-4 m/s 

90 km 

1000 m 

Latitude 

Ac 

Lx 

Lz 

28°S 

+7 m/s 

250 km 

1500 m 

33.5°S 

+ 10 m/s 

180 km 

400 m 

35.5°S 

+7 m/s 

130 km 

1200 m 

Latitude 

Ac 



28°S 

Om/s 



33.5°S 

-5 m/s 



35.5°S 

-2 m/s 




Core 

Depth 

1400 m 
1200 m 
1000 m 
1200 m 
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